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O ■ Abstract 

We study the Dirichlet boundary value problem for viscoelastic diffusion in polymers. We 
show that its weak solutions generate a dissipative semiflow. We construct the minimal trajec- 
tory attractor and the global attractor for this problem. 

^ ; 1 Introduction 

^ ■ 

c3 . The concentration behaviour for diffusion of penetrant liquids in polymers cannot always be described 
S by the Fickian diffusion equation 

Ou 

— = div{D{u)Vu), (1.1) 

^ ■ where u = u{t, x) is the concentration, which depends on time t and the spatial point x, and D(u) is 
O^l . the diffusion coefficient. The phenomena running counter to (1.1) include case II diffusion, sorption 
overshoot, literal skinning, trapping skinning and desorption overshoot [8l [TO j [TT | [T2 | [23 | l24" t 126] , 130]- 
There is a number of approaches which explain these non-Fickian properties of polymeric diffusion. 
They have much in common: they are usually based on taking into account the viscoelastic nature of 
polymers (cf. [17] and references therein) and on the possibility of glass-rubber phase transition (see 
e.g. [30J with some review). We are going to study the model which is due to Cohen et al. [7J |8| [T3]. 
• i— i ! The Fickian diffusion equation is replaced by the system 

' Ou 

— = DAu + EAa, (1.2) 
at 

da du 

— +p(u, a)a = fxu + u— . (1.3) 

Here the second variable a(t,x) is introduced (it is called stress), D and E are the diffusion and 
stress-diffusion coefficients, resp., /i and v are non- negative constants, and the scalar function (3 is 
the inverse of the relaxation time, for instance, (3 can be [8] taken in the following form: 

(3 = (3(u) = \{Pr + Pg) + \{Pr- Pa) tanh(^^) (1.4) 



ON 

in 
d 

o> 
o 



where p R ,p G ,5,u RG are positive constants, p R > Po% 



1 Formula (1.4) describes the following peculiarities of the processes under consideration. The polymer network 
in the glassy state (low concentration area) is severely entangled, so (3 is approximately equal to some small flc- In 
the high concentration areas the system is in the rubbery state: the network disentangles, so the relaxation time is 
small, and its inverse is close to (3r > @q- The glass-rubber phase transition occurs near a certain concentration urg- 
However, we assume that (3 also depends on stress, cf. ffi [Til 126] . 
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Well-posedness issues for initial-boundary value problems for systems of viscoelastic diffusion 
equations have been studied in (TJ (2J [15l [251 [26l [28] , see [261 EZ] for brief reviews. These results include 
investigation of system (1.2), (1.3) and more general settings (diffusion with variable coefficients). Let 
us only recall the main results on global (in time) solvability: strong solutions exist globally for v — 
and D = E in the one-dimensional case pQ, and for suitable non-constant stress-diffusion coefficient, 
but not for all initial and boundary data |15j : global existence of weak solutions for the Dirichlet and 
Neumann problems in the general setting with variable coefficients in the multidimensional case is 
shown in [26] and |28j . resp., without restrictions on the initial and boundary data. 

Let us also mention here a study of a system obtained from (1.2)— (1.3) by some simplification, 
in [19] . and paper p3], which touches upon some long-time behaviour issues for a free boundary 
problem for a polymeric diffusion model based on Fick's law. A result on long-time behavior (not in 
the "attractor framework") of the general second boundary value problem can be found in [28] . 

In this work we are interested in the long-time behaviour of the solutions to Dirichlet initial- 
boundary value problem for system (1.2), (1.3). We show that the weak solutions generate a semiflow 
on a suitable phase space with L 2 -topology (this means that there is a unique solution for any data 
from the phase space, and the solution semigroup is continuous in t and x). However, it is not 
clear whether this semigroup is asymptotically compact, and the phase space is not complete, so this 
impedes proving of existence of the usual global attractor for this semigroup in this phase space. A 
possible way out is to use the concept of minimal trajectory attractor. Thus, we construct a minimal 
trajectory attractor, which generates some generalized global attractor for the weak solutions of the 
problem in the (completed) phase space. 

The theory of trajectory attractors was created by G. Sell, M. Vishik and V. Chepyzhov [H 
El [20], in order to construct an attractor to weak solutions of the 3D Navier-Stokes equation. A 
generalization of this approach with a related notion of minimal trajectory attractor may be found 
in [311 [29] : it is applicable when the system lacks continuity properties or invariance of the trajectory 
space with respect to time shifts. In both theories, the trajectory attractor generates some generalized 
global attractor in the phase space. This global attractor has many usual properties of attractors, 
but its invariance may be shown only under additional conditions (see [31} Section 4.2.7]). The 
idea of trajectory attractor was slightly criticized by Sir J. Ball [I], for the evolution of the original 
system is not explicitly involved in its definition. However, an example [3Tj Remark 4.2.13] shows 
that the minimal trajectory attractor (and the corresponding global attractor) can well characterize 
the long-time behaviour of the system, even when the usual global attractor does not exist. That 
example also illustrates that such situations may appear, in particular, for problems with uniqueness 
of solutions (i.e. when there exists a solution semigroup). It is a rather unexpected fact because the 
original theory of trajectory attractors was developed for the problems where the uniqueness is not 
proved or is absentH. Similarly, in this paper we apply the theory of minimal trajectory attractors to 
a problem with uniqueness. However, in our case, we cannot insist that there is no attractor of the 
semigroup (semiflow), since this is unknown. 

Our paper is organized in the following way. In Section 2, we introduce the required function 
spaces. In Section 3, we give a weak formulation of the initial-boundary value problem for system 
(1.2)-(1.3) with existence, uniqueness and regularity results (Theorem 3.2, Remark 3.4). In Section 
4, we recall the basic issues of the classical and minimal trajectory attractor theories, and construct a 
dissipative semiflow generated by weak solutions of (1.2)-(1.3) (Theorem 4.10). In Section 5, we show 
that the trajectory space generated by the weak solutions possesses a minimal trajectory attractor 
and a global attractor (Theorem 5.2). 

2 Trajectory attractors for problems with uniqueness were investigated in [5] only as an intermediate step on the 
way to usual global attractors of semigroups. 
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2 Function spaces and related notations 

Lp(fi), W™{Sl), H m (n) = W^(Q) (m G Z,l < p < oo), # m ( fi ) =W" ( m e N ) are ' as 

usual, Lebesgue and Sobolev spaces of functions defined on a bounded open set (domain) Q C M n , 
fi£N. The scalar product and the Euclidian norm in L 2 {Vt) k = L 2 (n,R k ) are denoted by (it, v) and 
\\u\\, respectively (k is equal to 1 or n). In Hq(Q), we use the following scalar product and norm: 
(u,v)i = (Vtt, Vf), || u ||i = ||Vu||. We recall Friedrichs' inequality 

||u|| < #n||u||i- (2-1) 

Let ^(n) denote the topological subspace of L 2 (ST) consisting of functions from f/^Q). 
The Laplace operator A : Hq(Q) — > iJ _1 (fi) is an isomorphism. Therefore 

A -1 : if -1 (f2) — > i?o(^) (2.2) 

is also an isomorphism. Set X = X(Q) = A^^H^Q)). The scalar product and norm in X are 
{u,v) x = (A«,Ar))i, ||ii||x = ||Au||i. 

As usual, we identify the space i/ _1 (f2) with the space of linear continuous functionals on Hq(Q) 
(the dual space). The value of a functional from i7 _1 (f2) on an element from Hq(£1) is denoted by 
(•, ■) (the "bra-ket" notation). The scalar product and norm in i7 _1 (f2) are (it, = (A^u, A -1 d)i, 
||w||-i = ||A _1 m||i. Note that 

(u, Au)_i = u G iT -1 ^)^ e flj(fi). (2.3) 

The symbols C(J;E), L 2 {J\E) etc. denote the spaces of continuous, quadratically integrable 
etc. functions on an interval J C M. with values in a Banach space E. 

Let us remind that a pre-norm in the Frechet space C([0, +00); i?) may be defined by the formula 

+00 

IM|c([0,+oo);£) = ZJ 2_ Y 



M|C([0,t];g) 
( _j -I" |M|c([0,i];E) 



If E is a function space (L 2 (fl), H m (fl) etc.), then we identify the elements of C(J; E), L 2 (J~; E) 
etc. with scalar functions defined on J x Q according to the formula 

u(t)(x) = u(t,x), t G J, x G Q. 

We shall also use the function space (T is a positive number): 

W = W(Q, T) = {ue L 2 (0, T; H^(0)), u G L 2 (0, T; tf" 1 ^))}, 

IMIw = \\ u \\l 2 {0,T-HI{Q,)) + Il w lUaCO.Tii?- 1 ^)); 

[3~H Corollary 2.2.3] implies continuous embedding W C C([0, T]; L 2 (f2)). Moreover, 

<«',«> = ~|H| a , t;G W. (2.4) 

We also denote 

W /oc (fi,+oo) = {uG C([0,+oc);L 2 (Q)), W | [0)T ] G W(fl, T) VT > 0}, 

F/ oc (0, +00; Hl{tt)) = {ue C([0, +00); ^(fi)), M | [0 , T] G ^(0, T; ^(fi)) VT > 0}. 

We use the notation | • | for the absolute value of a number and for the Euclidean norm in W 1 . 
The symbols C stands for various positive constants. 
Note that 

(u,v) = (u,v), u G L 2 (Q),v G (2.5) 
and, since if -1 (f2) C L 2 (Q), one has 

H_i<C||m||. (2.6) 
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3 Basic properties of the boundary value problem 

We study the diffusion of a penetrant in a polymer filling a bounded domain C R n , n e N, which 
is described by the following boundary value problem: 

— = DAu + EAa, (t, x) G [0, oo) x tt, (3.1) 
Oct du 

— + P (u, o)a = fm + v—, (t, x) G [0, oo) x fi, (3.2) 
at at 

u(t, x) = <p(x), (t, x) G [0, oo) x dQ. (3.3) 

Here u = u(t,x) : [0, oo) x Q — ► R is the unknown concentration of the penetrant (at the spatial 
point x at the moment of time t), a = a(t, x) : [0, oo) x Q — > R is the unknown stress, <p : dfl — > R is 
a given boundary condition, fx, u, D and E are positive constants^, (3q : R 2 — > R is a given function, 

> Po(u, a) > 0g > 0; Q C R n is supposed to be a bounded open set such that X(Q) C Ty^fi) 
for some po > 2 (cf. [26]). For definiteness, we assume that the boundary condition for the stress is 
also prescribed: 

a(t, x) = 0(x), (t, x) G [0, oo) x dQ. (3.4) 
Equation (3.2) yields the following relation between the boundary conditions <fi and ip: 

P (<p, 4>)<j) = (J>¥>, x G dQ. (3.5) 
W.l.o.g. the functions ip and i/j are defined on Q. 

We assume that the functions j3 , 4> and -0 are C 2 — smooth. Moreover, for simplicity, let /3 (u, a) = 
Poo > 0g f° r large \u\ + \a\ (this assumption is admissible in the considered model, see Remark 3.6 
below). 

Set 

v(t, x) = u(t, x) — <p(x), w{t, x) = a(t, x) — 4>(x), 
(3{x, v, w) = (3 (v + ip(x),za + <f)(x)), 
h(x) = DAtp(x) + EA(f)(x), 
g(x, v, w) = ^>p{x) - p (v + ip(x),w + (f)(x))(p(x). 
Then we can rewrite (3.1)- (3.4) in the following form: 

dv 

— = DAv + EAw + h, (3.6) 
at 

dzu . dv 

— + p{x, v,w)zu = g(x, v, w) + fiv + u—, (3.7) 

v\dn = Mdn = 0. (3.8) 

This form of the studied problem will be useful below. However, in order to set the problem 
finally, we need the variable r(t, x) = zu(t, x) — pv(t, x). We denote 

d = D + uE, 

7 (x, v, r) = fiv — P(x, v, r + vv)t — is/3(x, v, r + vv)v + g(x, v , r + vv). 



3 The case fi = ("the Maxwell model" 0[3O]) is admissible as well. 
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Then problem (3.6)-(3.8) becomes 

dv 

— = dAv + EAt + h, (3.9) 
dr 

-^ = l{x,v,r), (3.10) 

v\an = r\ dn = 0. (3.11) 
It can be completed with the initial condition: 

v(0, x) = v (x), r(0, x) = t (x), x Eil. (3.12) 

Definition 3.1. A pair of functions (v,t) from the class 

v G W loc (Q, +oo), t G Hl oc (0, +oo; H^Q)) (3.13) 

is a weak solution to problem (3.9)-(3.11) if equality (3.9) holds in the space i7 _1 (f2) for a. a. t G 
(0,+oo), and (3.10) holds in H 1 ^) a.e. on (0,+oo). 

Theorem 3.1. Given vq G L 2 (Q) an d T o e Hq(Q), there exists a unique weak solution to problem 
(3.9)-(3.11) which belongs to (3.13) and satisfies (3.12). 

Note that (3.12) makes sense due to the embeddings W(fl,T) C C([0,T]; L 2 (fi)), 
H^O, T; Hq(Q)) C C([0, T]; H^(Q)), T > 0. 

Proof. Let us calculate the gradient of 7(2;, t>, r): 

V^(x, v, r) 

dp 

= liVv — P(x, v, t + vv)Vt — — (x, v, t + vv)t 

ox 

dp , \ „ 

— — {X, V, T + VVjTWV 

ov 

dp , ^ „ dp, . „ 

— v— — (x, v , r + uvjtvv — — — (x, r + vvjtvt 
Ovj Ovj 

nt dP, . 

— l'P{X, V, T + VV)\/V — U—[X, V, T + VV)V 

OX 

dP , x 2 ^ / x . . _ 

— IS-— (X, V, T + VV)V\V — V — — (X, V, T + VV)V\V — V— — (X, V, T + VV)V\T 

OV ovj ovj 

dg . . dg . dg , dg . 

+ — fx, f , r + z/f ) + — (x, f, r + uvjVv + ^- — fx, v, r + uvWv + — — fx, v,r + uvWt. 

OX OV OVJ ' OVJ 

Compactness of fl gives boundedness of 4> and tp. Hence, the function P(x,v,zu) is equal to P^ 
for large \v\ + \vj\. This implies that the gradient V7(x,t>,r) can be expressed in the form 

hi(x, v, t)Vv + h 2 (x, v, r)Vr + 9(x, v, r), 

where 

\hi(x,v,r)\ + \h 2 {x,v,r)\ < K h , 
\6{x,v,r)\ < K e (\v\ + |r| + 1) 

with some constants Kh, Kg. 
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Therefore, by Theorem 3.1], there is a pair of functions (v, r) from (3.13) which satisfies (3.9), 
(3.12) and 

A^ = A[ 7 (x,v,t)], (3.14) 

in ii/" _1 (f2) a.e. on (0,+oo) (more precisely, that theorem gives existence of solutions on finite time 
intervals, but then the solution on the positive semi-axis can be constructed step by step in a standard 
way, see 



Let us show that (v, r) satisfies (3.10). Since % G H£(0) for t > 0, it suffices to prove that 
j(x, v, t) G Hq{VL) for a.a. t > 0. Observe that j(x, 0, 0) = due to (3.5). The above representations 
of 7(x, u, r) and Vj(x,v,t) imply j(x,v(t,x),T(t,x)) G H 1 ^) (for a.a. t > 0). Let v m ,r m be 
sequences of smooth functions with compact supports in Q, v m — > v(t),r m — > r(t) in if 1 (f2). Then 
j{-,v m ,T m ) G Hq(Q). It remains to observe that j{-,v m ,T m ) — > 7(-, u(t), r(£)) weakly in if 1 (f2). 
Really, by Krasnoselskii's theorem [IB"! [22] on continuity of Nemytskii operators we have 

l(-,v m ,T m ) -> 7(-,v(t),r(t)),/i 1 (-,w m ,r m ) -> /h(-, u(t),r(t)), 

h 2 (-,v m ,T m ) -> /i 2 (-,i;(t),r(t)),6)(-,t; m ,T m ) -> 0(-,v(t),r(t)) 

strongly in L 2 (fi). Moreover, due to boundedness of /ii and /i 2 , w.l.o.g. we may assume that 
) -> /n(-, V (t),r(t)), /i 2 ( 

'^mi' T ra) > h 2 [- , v (t) , t (t)) *-weakly in Lcaiyi). But Vl> m > 
Vf (t), Vr m —> Vr(t) strongly in L 2 (f2). Therefore, 

hi(-,v m ,r m )Vv m -> /ii(-,v(t),r(t))Vw(t),/i 2 (-,f m ,r m )Vr m -> /i 2 (-, r(i)) Vr(t) 

weakly in L 2 (f2). 

Remark 3.1. It seems that a more profound reasoning of this kind proves that for any Nemytskii 
operator M : H l {tt) -> tf 1 ^) one has Af(£) - JV(0) G ^(fl) for all £ G #o(^)- 

It remains to prove uniqueness. If 07 = r+z/w, then tu G IV/ oc (fi, +00), and the pair (v, w) satisfies 
(3.6), (3.7). It suffices to show uniqueness of the pair (v,w). Let (vi,wi), (1)2,^2) be solutions of 
(3.6), (3.7) in the class Wi oc (fl, +00) x Wi oc (fl, +00) with the same initial conditions. Let us denote 
w = Vi — v 2 , £ = vo\ — W2- Then 

w' = DAw + EA£, (3.15) 

£' + (3(x, v 1 ,wi)zui - (3{x, v 2 , vj 2 )w 2 = g(x, v x , wx) - g(x, v 2 , w 2 ) + fiw + vw' . (3.16) 

Calculate the i/ _1 (f2)-scalar product of (3.15) and \m + vw' for a.a. t G (0, 00), and take (2.3) 
and (2.5) into account: 

fj,(w', w)-i + u(w', w')-x = —/jbD(w, w) — (J>E(£, w) — vD (w' , w) — vE (it/, £) . (3-17) 

Take the "bra-ket" of (3.16) and E£ for a.a. t G (0,oo), and remember (2.5): 

E (£', £) + E(fi(x, vi, zui)mi - 0(x, v 2 , w 2 )w 2 , £) = 

E(g(x, vx, wi) - g(x, v 2 , w 2 ),() + fiE(w, f ) + vE (w', f ) . (3.18) 
Adding (3.17) and (3.18), and omitting the second term, which is positive, we conclude: 

MINI 2 -! m ||2 vDd\\wf Ed\\jf 
2 dt ^ 11 " j " 2 dt 2 dt 

< E((3 (v 2 + v?,ro 2 + 4>){w 2 + 0) - + m + 4>){wi + (/>),£) 

<C(|K|| + |h||)||e||. (3.19) 
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The last inequality follows from boundedness of dlffoKo-) 0- ] an( j 9\M^e\^ Lagrange's theorem and the 
Cauchy-Buniakowski inequality. 
Integration from to t yields 

f 



£HI 2 -i + A^ / IK fl )|| 2 ds + ^|Hl a + flieir 





t 



<C l(U(s)\\ 2 +\Ms)r)ds. 

Hence, 

\Ht)\\ 2 + \mW <(' ll\\V«)f + \\<rl«)f)<l«. t > U. (.3.20) 



Thus, by the Gronwall lemma, £ = w = 0. □ 

Remark 3.2. Weak solutions (v(t),r(t)) to (3.9)-(3.11) belong to H^(tt) 2 for all t > 0. It follows 
from a simple regularity result for react ion- diffusion equations. Consider the problem 

dv 

— -aAv = /, (3.21) 

v\ d n = 0, (3.22) 

v(0,x) = v (x), xeVt. (3.23) 

Given a > 0, v G L 2 (fi) and / G 1/2(0, T; i7 _1 (fi)), there exists a unique weak solution v G T), 
T > 0. Since W{Q,T) C C([0, T], L 2 (fi)), the operator 

T : L 2 (fi) x L 2 (U } T;H-\ti)) -> L 2 (fi), T(uo,/) = u(*),0 < t < T, 

is well-defined. Since r G H x (0, T\ Hq(Q)), it suffices to apply the following lemma to equation (3.9). 

Lemma 3.1. The operator T transforms L 2 (£l) x H 1 ^, T; if _1 (f2)) into ^(fi). 

Proof. W.l.o.g. a = 1. The solution v can be considered as the sum V\ + v 2 of the weak solutions to 
the following problems: 

^ - At* = /, (3.24) 

ui|an = 0, (3.25) 

t; 1 (0) = -A- 1 /(0), (3.26) 

^ - At; 2 = 0, (3.27) 

v 2 |an = 0, (3.28) 

u 2 (0) = Uo + A- 1 /(0). (3.29) 

i 

It is easy to see that V\(t) = — A _1 /(0) + J v 3 (s) ds, where v 3 G W(Q,T) is determined by the 

o 

problem 

v' 3 - Av 3 = /', (3.30) 

v 3 \ 9 n = 0, (3.31) 

v 3 (0) = 0. (3.32) 

Obviously, v 3 G L 2 (0,T; H^(Q)), so ^ G H\0,T;H^n)) C C([0,T]; -Hj(fi)). Moreover, see e.g. 

P Proposition XV.3.5], u 2 (*) G ^(fl). □ 
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Remark 3.3. There is no essential loss in generality of the model in assuming that @o(u, a) is equal 
to some constant (3^ for large \u\ + |er|. In fact, physically, \u(t, x)\ < 1 (u is the concentration, so it 
cannot exceed 100%). Let q = a — vu. Then (3.2) yields (cf. [26l [27] ) 

q(t,x) = ^(0,x)exp ( - J x), vu(£, x) + q(£, x)) d£ 



Hence, 



q 



+ j exp I J p (u(€,x),vu(£,x) + q(q~,x))dq~ 
o \t / 
x [ft — v(3 (u(s, x), vu{s, x) + q(s, x))] u(s, x) ds. 



t 

(t, x) I < e- tl5c \q(0, x)\ + + u/3 R ) [ e (s "' )/3G ds < e~^ G \q(0, x) \ + ^ + V ^ R . (3.33) 







0, 



G 



Thus, if |^(0, x) I is uniformly bounded, q is also bounded: 

\q{t,x)\ < C,Vt > 0,x e O. (3.34) 

Moreover, (3.33) implies that long-time behaviour of q(t,x) is bounded for any \q(0, x)\. Hence, 
q (and, therefore, a) is bounded for typical regimes which may be observed in reality. Thus, the 
relaxation time (and its inverse (3q) can be experimentally determined^ only for bounded u and a, 
whereas "at infinity" we can choose it at discretion, for instance, we can let (3q{u, a) = 0^ for large 
\u\ + \a\. 



4 Semigroups, semiflows, trajectory attractors and global 
attractors 

Let us recall some basics of the attractor theory. Let E be a metric space. 

Definition 4.1. A family of mappings St ■ E — > E, t > 0, is called a semigroup if So is the identity 
map / and 

S t oS s =S t+s (4.1) 

for any t, s > 0. 

Definition 4.2. A set P C E is called attracting (for S t ) if for any bounded set B C E and any 
open neighborhood W of P there exists h > such that StB C W for all t > h. 

Definition 4.3. A set P C E is called absorbing (for S t ) if for any bounded set B C E there is 
h > such that for all t > h one has S t B C P. 

Definition 4.4. A set A C E is called invariant (for St) if 

S t A = A 

for any t > 0. 

4 e.g. in form (1.4). 
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Definition 4.5. A set A C E is called a global attractor (of St) if 

i) A is compact; 

ii) A is invariant for iS t ; 

iii) A is attracting for S t . 

If there exists a global attractor of S t , then it is unique (see e.g. [2U Corollary 4.1.1]). 

Definition 4.6. A semigroup St : E —>■ E is called dissipative if there is a bounded absorbing set. 

Definition 4.7. A semigroup St '■ E — > E is called asymptotically compact if, for any bounded 
sequence y m <E E and any sequence of numbers t m — > +oo, the sequence <S fm (y m ) contains a converging 
subsequence. 

Definition 4.8. A semigroup «S t : E —>■ E is called a semiflow if the map 

(t,y) i-> <S t (?/) 

is continuous from [0, +oo) x E to E. 

The next result follows e.g. from [H Theorem 3.3, Corollary 4.3]: 

Theorem 4.1. A semiflow St : E —>■ E has a global attractor A C E if and only if it is dissipative 
and asymptotically compact. If E is connected, then A is a connected set. 

In order to describe the dynamics of weak solutions for problem (3.9)-(3.11), one may put E = 
L 2 (£l) x L\{VL), and define the semigroup St '■ E — > E in the standard way: if y = (v , r ), and (v, r) is 
the corresponding weak solution of (3.9)-(3.12), then we set S t (y) = {v(t), r(t)). Since weak solutions 
belong to C([0,T]; L 2 (fl)) x C([0,T]; if 1 (fi)) for all T > 0, the map t t— > S t (y) is continuous for each 
y G E. On the other hand, a reasoning similar to the proof of uniqueness in Theorem 3.2 shows that 
the map S t : E — > E is continuous uniformly with respect to t G [0,T] for all T > 0. Then the map 
(t,y) i — ► St(y) is continuous, and, taking into account Lemma 5.1 (see below), we arrive at 

Theorem 4.2. St ■ E —>■ E is a dissipative semiflow. 

However, it seems to be very hard or impossible to establish asymptotic compactness of this 
semiflow since the space E is not complete (and this is usually important, cf. [H]). If we try to 
change the phase space E and to find so-called (E, F)-attractors [SI EI], i-e. attractors which attract 
in the topology of some space F, which is weaker than the one of E, then we need continuity of S t in 
this weaker topology, which is not clear. Moreover, it seems (due to the properties of the generalized 
attractor which we construct below) that the attractor can contain non-differentiable elements of 
L 2 (f2), and then Remark 3.4 implies that it cannot be invariant. Thus, we are going to use the 
theory of trajectory attractors, so let us briefly describe the required notions. 

Let E and E be Banach spaces, E C E , E is reflexive. Fix some set 

H + C C([0, +oo); E ) n L^O, +oo; E) 

of solutions (strong, weak, etc.) for any given autonomous differential equation or boundary value 
problem. Hereafter, the set 7i + will be called the trajectory space and its elements will be called 
trajectories. Generally speaking, the nature of 7i + may be different from the just described one. 

Definition 4.9. A set P C C([0, +oo); E ) D £oo(0, +oo; E) is called attracting (for the trajectory 
space H + ) if for any set B C 7i + which is bounded in Loo(0, +oo; E), one has 

SUp inf \\T(h)u - V \\c([0,+oa);E ) 0. 
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Here T(h) stands for the translation (shift) operators, 



T(h)(u)(t) = u{t + h). 
Note that T{h) is a semiflow on C([0, +00); E ). 

Definition 4.10. A set P C C([0, +00); Eq) R Loo(0, +00; £?) is called absorbing (for the trajectory 
space H + ) if for any set B C which is bounded in Loo(0, +00; E), there is h > such that for all 
t > /i: 

T(t) J B C P. 

Definition 4.11. A set U C C([0, +00); £ , )nL oo (0, +00; i?) is called the minimal trajectory attractor 
(for the trajectory space Ti + ) if 

i) U is compact in C([0, +00); E ) and bounded in Loo(0, +00; E 1 ); 

ii) T(i)W = U for any t > 0; 

iii) U is attracting in the sense of Definition 4.11; 

iv) U is contained in any other set satisfying conditions i), ii), iii). 

Definition 4.12. A set A C E is called the global attractor (in Eq) for the trajectory space 7i + if 

i) A is compact in Eq and bounded in E; 

ii) for any bounded in Loo(0, +00; E) set B C 7i + the attraction property is fulfilled: 

sup inf ||tt(t) — v\\ Eo — > 
ueB 

iii) ^4 is the minimal set satisfying conditions i) and ii) (that is, A is contained in every set 
satisfying conditions i) and ii)). 

Theorem 4.3. (see Wl\ Corollary J^.2.1, Lemma J^.2.9]) Assume that there exists an absorbing 
set P for the trajectory space 7i + , which is relatively compact in C([0, +00); Eq) and bounded in 
-^oo(0, +00; E). Then there exists a minimal trajectory attractor U for the trajectory space 7i + . 

The structure of minimal trajectory attractors is discussed in [2H Chapter 4] (see also [B])- In 
particular, the minimal trajectory attractor contains the set of those solutions to the considered 
problem that can be continued to the whole real axis being uniformly bounded in E and continuous 
with values in E Q ; however, some extra elements may appear. The structure of global attractors is 
determined by the structure of minimal trajectory attractors: 

Theorem 4.4. (see I31\ Theorem J^.2.2}) If there exists a minimal trajectory attractor 14 for the 
trajectory space 7i + , then there is a global attractor A for the trajectory space 7i + , and for allt > 
one has A = {y(t)\y G U}. 

5 Attractors for the polymeric diffusion problem 

We are going to construct the minimal trajectory attractor and the global attractor for problem 
(3.9)-(3.11). Let us choose L 2 (fi) 2 as the space E and the space H S (Q) 2 as the space Eq, where 
5 G (0, 1] is a fixed number. The trajectory space 7i + is the set of all weak solutions to (3.9)-(3.11) 
from class (3.13). It is contained in L^O, +00; E) due to Lemma 5.1 (below) and in C([0, +00); E ) 
(clearly). 

Dissipativity of problem (3.9)-(3.11) is stated by the following result: 
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Lemma 5.1. Let (v,r) be a weak solution to (3.9)-(3.11). Then 

^IK^f + flkW + ^f + flKOf-x^ 



e- 7i lk(0)|| 2 + -||r(0) + ^(0)|| 2 + 5||t;(0)||M + V (5.1) 
/or a// 1 > 0, where T and 7 are some fixed positive numbers, which are independent of v, r and t. 

Proof. The pair (v,w = r + vv) satisfies (3.6), (3.7). Take the if _1 (f2)-scalar product of (3.6) and 
\iv + vv' and the "bra-ket" of (3.7) and Ew for a. a. t G (0, 00), and add the results (cf. proof of 
Theorem 3.2): 

- " ," + v \\v 2 i + W^ M H V - H ^-r- - 

2 dt 11 11-1 p 11 11 2 dt 2 dt 

+E(P(x, v, w)w, w) = (h, fiv + vv')_i + E(g(x, v, w),vj). (5.2) 



Hence, 



//d|H| 2 ! 2 2 vDd\\v\\ 2 Ed\\w\\ 2 



+ HKII-i + nD\\v\\* + — — + 



2 dt " 11-1 ^ " " 2 dt 2 dt 

+£/3 G |M| 2 < C(|M|-i + IKH-i + |M|) < C(\\v\\ + ||v'||_i + H|). (5.3) 

The Cauchy inequality for scalars can be written in the form Cr] < erf + Thus, (5.3) implies 
A* AH-i , /*Q„ „ a , yDd|Ml a , ^ d|| CT |l 2 < E/3 G|| 

; \\v\\ H ; ; W < O. 5.4 

2 dt 2 11 11 2 dt 2 dt 2 11 11 - v ; 

Let x (t) = f ||«<*)ll— 1 + f IKOII 2 + f H(f)|| 2 . Obviously, 

^II^)H 2 + ^II^)I| 2 >7XW 

for some 7 > 0. Thus, ^(t) + 7X(*) < c , so x(*) < e~ 7 *x(0) + (1 + 7~ 1 )C by Lemma II.1.3], and 
(5.1) follows. □ 

The main result of this section is 

Theorem 5.1. The trajectory space 7i + possesses a minimal trajectory attractor and a global attrac- 
tor. 

Proof. Due to Theorems 4.15 and 4.16, it suffices to find an absorbing set for the trajectory space 
7i + , which is relatively compact in C([0, +00); E ) and bounded in Loo(0, +00; E). Consider the set 
P of weak solutions to (3.9)-(3. 11) which satisfy the estimate 

^-\Ht)\\ 2 + §\\r(t) + MtW + %Mt)\\-i < 2r >Vt > 0. (5.5) 

It is an absorbing set for the trajectory space H + and is bounded in L oo (0,+cx3; E). By (3.9), the 
set {v', (v, t) G P} is bounded in L^O, +00; H~ 2 (tt)). Moreover, (3.10) yields that {r', (v, r) G P} 
is bounded in Loo(0, +00; L 2 (f2)). The embedding E C E is compact. By [211 Corollary 4], the set 
{y\[o,M],y G P} is relatively compact in C([0, M];E ) for any M > 0. This implies (cf. |3l, p. 183]) 
that P is relatively compact in C([0, +00); i? )- □ 
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